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We consider a network model of snake states to study the localization problem of non-interacting 
fermions in a random magnetic field with zero average. After averaging over the randomness, the 
network of snake states is mapped onto M coupled SU(2iV) spin chains in the N — * limit. The 
number of snake states near the zero-field contour, M, is an even integer. In the large conductance 

2 

limit g = M^-j- (M S> 2), it turns out that this system is equivalent to a particular representation 
of the U(27V)/U(iV) x U(iV) sigma model (iV — » 0) without a topological term. The beta function 
/3(1/M) of this sigma model in the 1/M expansion is consistent with the previously known /3(g) of 

l/") ' the unitary ensemble. These results and further plausible arguments support the conclusion that 

ON , all the states are localized. 
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There has been much recent interest in the quantum motion of a charged particle in a random perpendicular 
magnetic field in two dimensions. Theoretically, this problem arises from the study of gauge theories of strongly- 
correlated electronic systems. An example is the Chern-Simons theory for the half-filled Landau level In a 
mean-field treatment, one obtains an effective theory of free fermions moving in a weak effective magnetic field due 
to the overall cancellation of the external field by the Chern-Simons field. This effective field has zero average but 
contains a random static component induced by the inhomogeneous electron density. Another example is the gauge 
theory of the doped Mott insulators J3|,[| . In this case, slow fluctuations of the gauge field can be approximated as a 
\0 ! static random field. 

Single-particle motion in a random magnetic field with zero average is also interesting in the context of localization 
theory. One might expect that a system with random magnetic fields should belong to the same universality class as 
other disordered systems with broken time-reversal symmetry and zero Hall conductance. This is the unitary class 
' for which all states are localized, according to conventional scaling theory ||. In field-theoretic terms, the system 
ON ■ would be described at large distances by a non- linear sigma model (without a topological term). However, Zhang and 
Arovas Q have recently suggested that long-range interactions between topological densities in the non-linear sigma 
model, induced by local fluctuations of <r xy , may give rise to a novel derealization transition when the conductance 
Q ■ of the system at short distances exceeds a critical value of the order of e 2 /h. 

In this paper, we will present an analytic approach to this random flux problem in the limit of smooth disorder, 
based on the random network model developed in Ref. Q| and the method introduced by D.H. Lee and coworkers 
P-|TT| for the discussion of the Chalker-Coddington network model in the context of the quantum Hall transitions. 
We find that the transport properties of the random-magnetic-field problem can be described in the large- conductance 
limit by a non-linear sigma model without a topological term. Therefore, we argue that the localization properties of 
. i | this system are the same as that of the conventional unitary ensemble, suggesting that all states are localized. Similar 
conclusions have also been reached for the case of weak and short-rang ed disorder HQ. 

There have been several numerical studies of this problem with conflicting conclusions. The difficulty in the 
interpretation of the numerical data arises from the fact that the localization length increases rapidly as a function 
of energy as one approaches the band center. Some authors |]l5|-|l8| have argued that there exists a range of energies 
for which the eigenstates are extended, while others |7p|,|l9[ suggest that all the states are localized. Sugiyama and 
Nagaosa (l9) have observed that the numerical data are consistent with the conventional single-parameter scaling of 
the conductance in the unitary universality class. Although their numerical data are restricted to bare conductances 
below 0.2e 2 //i, we will argue that the conventional scaling function continues to apply in the large-conductance limit 
of this problem, in contrast with the Zhang- Arovas theory. 

We consider in this paper the case of a random magnetic field with zero average (B(r)) = and a non-zero variance 
(B 2 (r)) — Bq. We will work in the semiclassical limit in which the correlation length A of the random magnetic 
field is large compared to the typical magnetic length 1$ = (hc/eB^) 1 / 2 . In this limit of smooth disorder, an electron 
moves along contours of constant magnetic field in a direction determined by the sign of the field gradient on the 
contour. This chirality reflects the broken time-reversal symmetry of the problem. Since we are interested in the 
possibility of derealization, we can focus on contours which percolate geometrically across the system. These are the 
B = contours where the classical motion of the electron describe snake- like trajectories To simplify further, we 
represent these zero-field contours as links on a square network, carrying a fixed number, M, of the electronic "snake 
states" . The nodes of the network represent the saddle points of the magnetic field, where a particle is scattered to 
the right or left. This scattering, which is chaotic even in the classical limit, is modeled by random amplitudes chosen 
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from a distribution which ensures that the network as a whole has zero Hall conductance. Quantum interference due 
to these scattering processes is modeled by random phase shifts and the mixing of the snake states along the links. 

Details of the mapping of the semiclassical limit to the network model can be found in Ref. . We will emphasize 
here only the points of particular relevance to the present discussion. It should be noted that the electron states on 
the zero-field contours typically occur in pairs at a given energy. This can be seen most easily by considering the 
special case of a straight contour, across which the field switches abruptly between two values, ±£?o M]. Restricting 
attention to the energies between the lowest two Landau levels in the bulk, it can be shown that two modes propagate 
along the contour in the same direction, arising ultimately from symmetric and antisymmetric linear combinations 
of the Landau levels on either side of the contour. Similarly, at a higher energy, each bulk Landau level below that 
energy is associated with two states on the zero-field contour. Thus, each link of the network carries an even number, 
M, of channels propagating in the same direction. 

The conductance of the system at short distances is given by 

A schematic representation of the network is given in Fig. 1 for the two-channel case which was studied numerically 
in Ref. @. In this paper, we are interested in the limit of large conductance so that M 3> 2. It should be pointed 
out that the semiclassical picture of motion confined to contours of constant field breaks down when the energy of 
the particle is so high that the lateral extent of the wavefunction becomes comparable to the correlation length of the 
magnetic field. In terms of the number of channels, this condition for the validity of the semiclassical picture is given 
by M < M max = (E/E slopc ) 3 / 4 , where E slopc = -f^r(VA) 2/3 [§■ Nevertheless, we can obtain an arbitrary M max by 
choosing a sufficie ntly large correlation length. 

Following Ref. |10| , we represent the network model using the schematic picture depicted in Fig. 2 for the two- 
channel case. Let us take each zigzag line as a constant-x curve and measure the displacement along it by the 
coordinate y. We use the coordinate system where the saddle points are located at (x, y) = (n x ,n y ), where n x and 
riy are integers. Taking the continuum limit in the y-direction, the network can be represented by the Hamiltonian 

§m 

X j 

* i 

-^2 / d V^2 \tjk{x,y) i>j(x,y)ip k (x,y) +h.c. , (2) 

x J j<k 

where ipj (j — 1,---,M) is the fermion annihilation operator of the jth snake state. The scattering processes of 
the network are now modeled by tunneling (tj) between neighboring sites with the same channel index, as well as 
tunneling (tjj~) between different channels at the same x-site. We choose tj(x,y) and tjk(x,y) to be independent 
random complex numbers with Gaussian distributions. The random magnitudes of these matrix elements incorporate 
local fluctuations of Hall conductance explicitly. In addition, the random phases will give rise to interference effects 
which may lead to the localization of these states. The two-channel network studied in Refs. (tJJ^] correspond here to 
the case where |ii| and fa] are chosen such that there is no bias for a particle to be scattered to the right or left at 
each node (if the node is considered in isolation). Since the effects of changing the energy, E, can be absorbed into a 
dependence of the tunneling matrix elements, tj(x, y) and tjk(x, y), on E, we can set E = without loss of generality. 
Note that, as we change E, both of the tunneling matrix elements and the number of snake states, M, change as 
functions of E. 

In order to study transport properties, we have to consider the generating functional for the retarded and advanced 
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=E/<^EE 

x j P 

+ E / dy EE \pj( x >v) 1 PvA x + ij^pj^.y) + h-c] 

X J p 

+ ^2 J dy^2^2[t jk (x,y) -ip pj (x,y)ip pk (x,y) +h.c] . (3) 

Here, and V'pj are Grassmann variables. Also, p = +(— ) represents the retarded (advanced) sector, S p — sgn(p), 
and 77 is a positive infinitesimal number. 

Now we let tp P j — > ippj(iip P j) and — > —i'4>pj{'ippj) for even (odd) i's |lC(| . Then, the action becomes 

•5 = E / ^EE K -1 )* 7 ^ ^pi&yWpj&y) + 4>pj(x,y)dy'>Ppj(x,y)] 

x j p 

+E/ *EE N 1 '^ ^w^+^fe^fi+H 

X i p 

~*E / ^EE( _1 ) X [tjkfav) $ P1 {x,y)^pk{x,y) +h.c] . (4) 

X j</c p 

Next we replicate the action and average over tj and tjk to get 

5 = J^ / dy^2^2[{-l) x rfS p xp aj (x,y)ip aj (x,y) + $ a j(x,y)dyip aj (x,y)] 

x ^ j<fc a,/3 

where a = (p, n.) and n = 1, • • • , iV is the replica index so that a and /3 take 2iV values. 

Taking y as the imaginary time, we can see that the above action describes a many-body system of interacting 
fermions on M one-dimensional chains. We can further substitute the fermion operators with the generators of 
SU(2JV): 



S h = ^« " S «0 ™ E ^ 
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Thus, we obtain a theory of M coupled SU(2iV) spin chains, described by the Hamiltonian: 

X j 

+ EE J i ^[S,(x + l)S 3 {x)] 

X j 

+ EE JjkTriS^Skix)} , (6) 

X j<k 

where J, = (|t,| 2 ), Jjk = -{\tjk\ 2 ) and A is a 2iV x 2N diagonal matrix with A^ = 1 for a < N and A^ = -1 
for N < a < 2N. Also Tr [AS] = J2af3 f° r an y matrix A and B. Since iJ commutes with the local density 

rij(x) = ^aipajipaj of each chain j, Hilbert spaces corresponding to different {n,j(x)} decouple. Consider first the 
case of decoupled chains, described by the first two terms of the Hamiltonian of Eq. (0) |jTo| , p7|| . In this case, the 
ground state of each chain belongs to the Hilbert space where rij(x) = N for all x [E0[. In this Hilbert space, a 
particular representation for SU(2iV) is realized, characterized by a Young tableau with a single column of length N 
JTo| , p7t . In the above mapping from Eq. (||) to Eq. (||) , we drop constant terms expecting that the ground state of the 
Hamiltonian (j|) still lies in the Hilbert space where nj(x) — const, for all x. 

The above Hamiltonian gives nearest- neighbor antifcrromagnetic couplings (Jj > 0) in each spin chain j. On the 
other hand, for each x, the coupling between spins on different chains, j and k, is ferromagnetic (Jjk < 0). In terms 
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of the original problem, the localization corresponds to the existence of an energy gap in the spectrum of this spin 
chain (in the replica limit of TV — > 0). 

Before considering the limit of many coupled chains, it is worthwhile to consider first the system for SU(2) with 
two coupled chains only (M — 2 and N = 1). In this case, we have two antiferromagnetic spin-^ chains with 
ferromagnetic coupling (Ji 2 = K± < 0) between the chains (see Fig. 3). We can consider a larger parameter space 
where the antiferromagnetic coupling along each spin chain alternates from bond to bond between two values, K\ and 
K 2 . This means that each spin chain corresponds to a single-channel network model with non-zero Hall conductance. 
In order to observe the constraint of zero overall Hall conductance for the present problem, we require the alternation 
of the bond strengths on the two chains to be staggered. We are therefore coupling together two single-channel 
networks with Hall conductances of opposite signs. In other words, we have J\(x = even) = J 2 (x = odd) = K\ and 
J\{x — odd) = J 2 (x — even) = K 2 with Ki,K 2 > 0. We can now study the system as a function of the parameter 
A = (Ki — K<z)/(Ki + K 2 ). Since the question of localization is believed to be independent of microscopic details, it 
is interesting to ask whether the spectrum of the two coupled chains has an energy gap irrespective of the parameters 
A and K± (as long as the latter is finite). 

As already mentioned, setting K±_ = results in two independent single-channel networks. This has been studied in 
detail by D.H. Lee and coworkers ]9|-[iT|] in the context of the plateau transition in quantum Hall systems. Each of these 
networks should exhibit delocalization at A = 0. For finite |A|, the alternating coupling gives rise to dimerization, 
i.e., a spin-Peierls phase. Instead of changing A to destroy the gapless phase, we may also switch on the interchain 
ferromagnetic coupling. This is a relevant perturbation to the decoupled gapless spin chains. Therefore, in the case 
of A = 0, we expect the low-energy physics of the coupled system to be described by the opposite limit of strong 
coupling (I-Kjj ^ Ki — K 2 ) where one has a uniform antiferromagnetic spin-1 chain. It is known that this system is 
always in a Haldane phase p^j25[ |. 

It remains to consider the case of finite A where we couple two chains in spin-Peierls phases in such a way that 
their dimerizations are staggered. This case has been studied numerically with the conclusion that the coupling does 
not destroy the energy gap. Instead, the system is believed to cross over from the dimerized phase to a Haldane phase 
p6| as K± becomes much stronger than K\ and K 2 . Thus, we can see that there are no phase transitions along the 
boundaries of the A — K± plane. This suggests that there are also no phase boundaries within this A — K± plane. 
In other words, the existence of an energy gap in this system is indeed independent of these parameters of the mode, 
and it points to the absence of a delocalization transition in the network model. 

Let us return to the general case of arbitrary M and N. As in the SU(2) case considered above, the interchain 
ferromagnetic coupling between any two of M spin chains is a relevant perturbation in the renormalization-group 
sense. As the number M of spin chains becomes large, we expect the M SU(2A r ) spins at the same x-site to form 
a totally symmetric representation of SU(2iV), which corresponds to a Young tableau with M columns of length N. 
This strong-coupling limit can be regarded as a single SU(2iV) spin chain where the spin at each site is in this large- M 
representation. The residual coupling between the spins at different sites remain antiferromagnetic. 

In the large- M limit (as 77 — > 0) , this antiferromagnetic spin chain is equivalent to the following U(2iV)/U(A r ) xJJ(N) 
sigma model [pl JTlpTll 
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C = —Tr(d,Q) 2 + y^Tr [Qd^Qd v Q] , (7) 

where fj, = x,y, Q = m'Am, and u is a 2N x 2N unitary matrix. The second term is a topological term [p3| , [27f and 
satisfies 

J drdx e^Tr [Qd^QdvQ] = 16mm , (8) 

where m is an integer. Since M is always an even integer in our problem, the topological term is always irrelevant and 
can be dropped in the Lagrangian. Therefore, the system can be described by a U(2iV) /U(N) x U(iV) sigma model 
without a topological term. Zhang and Arovas j^] have argued that a term describing logarithmic interactions between 
topological densities may exist and lead to a gapless phase, similar to the Kosterlitz-Thouless phase. However, we 
have been unable to find any evidence for such behavior in this spin-chain mapping. In particular, for the N = 1 
problem of SU(2) spins, it is known that the spin chain has a gap in its spectrum for all even M (integer spins). 
Thus, one can see that we do not have the long-range interaction between topologocal densities in the case of N — 1. 
According to Zhang and Arovas [||, there should be a term describing the long-range interaction for arbitrary N. 
Therefore, unless the theory is highly nonanalytic in terms of N so that the replica limit is not well defined, we expect 
the absence of long-range interaction between topological densities even for arbitrary N. 

The beta function of the non-linear sigma model (without a topological term) in the N — > limit was evaluated by 
Hikami fl28|] and is given by 
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(3(1/M) ee d = 2(1/M) 2 + 6(1/Af) 4 + 0{(l/Mf) , (9) 

where L is the size of the system. Using the relation Eq. (Q) , we can obtain the beta function for the dimensionless 
conductance g = G/(e 2 /h): 

This is the known beta function for the dimensionless conductance in the unitary ensemble |l3|p8| ]. We therefore 
believe that all the states are localized and that the localization length in the large-conductance limit is given by 
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In summary, we have used a multi-channel network model to study the localization properties of non-interacting 
fermions in a smoothly varying random magnetic field with zero average. Using the replica trick, the network with 
an even number of channels is mapped onto an even number of coupled SU(2Y) spin chains in the N — > limit, 
where TV is the replica index. In the large-conductance limit where the number of channels is large, this system is 
equivalent to a particular representation of the U(2iV)/U(iV) x XJ(N) sigma model (N — > 0). Thus, we confirm that 
the beta function (3(1 /M) of this model is consistent with the known (3(g) of the unitary ensemble. This result, which 
suggests that all states are localized, agrees with perturbative calculations |l^ , |T4|]_f or random fields with short-ranged 
correlations, and should be contrasted with the prediction of Zhang and Arovas||6| that the system is critical for bare 
conductances greater than e 2 /h. 
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FIG. 1. A schematic illustration of a network in the problem of the random magnetic field. The two channels are denoted 
by solid and dashed lines. 

FIG. 2. A representation of the network in terms of chiral fermions with alternating directions in the case of M = 2. The 
two channels are denoted by solid and dashed lines. 

FIG. 3. A schematic representation of the two coupled spin chains (M = 2). Note that Jj(x) — Jj(x + 2) [Jj(x) = (\tj\ 2 )(x) 
and j = 1,2]. Also, Ji(x = even) = 3i(x — odd) = K\ and J\(x = odd) = 3i{x — even) = K2. J12 = (\tjk\ 2 ) = K± is 
independent of x. 
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